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A limit on the mass of the tau neutrino m^^ is derived from 4.5 x 10^ tau pairs produced in an integrated 
luminosity of 5.0 fb"^ of e^e" — > 7* t^t~ reactions at center of mass energies « 10.6 GeV. The measurement 
technique involves a two-dimensional extended likelihood analysis, including the dependence of the end-point 
population on m^^ , and allows for the first time an explicit background contribution. We use the decays r 5nvT 
and r ?i-k2-k'^Vt to obtain an upper limit of 30 MeV/c^ at 95% C.L., as well as a preliminary limit of 31 MeV/c^ 
from T Stttt'^Ut 
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1. Introduction 

We present two separate analyses of the mass 
of the i/r- The first analysis, as published in [Q, 
herein referred to as the 5h analysis, comprises a 
sample of the decays r —^ birvT and t 37r27r*^^'T-. 
These decay modes have a small overall branch- 
ing fraction but a significant probability of pop- 
ulating the hadronic mass end-point sensitive to 
m^. The second (preliminary) analysis, herein 
referred to as the Ah analysis, exploits the decay 
T — > Stttt'^Vt- Although the probability per event 
of being near the endpoint is much smaller in this 
sample, the much larger branching ratio makes it 
competitive. 

2. The CLEO Analysis Method 

The analyses use the two-dimensional hadronic 
energy versus hadronic mass spectrum, first used 
in , to obtain a limit on the mass of the . 
The selection criteria require well reconstructed 
tracks, as well as good quality tt^ candidates, 
when appropriate, recoiling against a well sepa- 
rated one prong tag. For the 5h analysis, the tag 
track is required to be an identified lepton, while 
in the Ah analysis, the tag must be consistent with 
either an electron, muon, tt or p. Stringent excess 
energy criteria reduce the possible backgrounds in 
both analyses. In the 5h analysis only a small re- 
gion near the endpoint is used for the fit. For that 



analysis, an extra Poisson coefficient is used in an 
extended likelihood, relating the number of events 
expected in the Fit Region at the endpoint as a 
function of neutrino mass to the number of events 
observed in data in a neutrino-mass-insensitive 
Control Region (see Fig0.) The fit region for the 
Ah analysis is chosen large enough to make the 
Poisson term unnecessary. For the first time in 
any Vr mass analysis , both the 5/i and Ah fits 
allow for an explicit background term to account 
for non-signal events. 

3. The 5h Sample 

The final event sample is summarized in Ta- 
bic |l|. The hadronic masses of the 5h samples 
are displayed in Fig |^. There are 207 events in 
the 37r27r" sample and 266 events in the Stt sam- 
ple. Backgrounds can come from either misre- 
constructed tau decays or from non-tau sources 
such as charm decays. The non-tau background 
estimator shape used in the fit is obtained from 
a data sample of events with tracks in one hemi- 
sphere with loose selection criteria and invariant 
mass above the tau mass recoiling against either 
57r or 3tt2tt^ candidates. Its normalization is pro- 
vided by the number of events in Fig || above the 
tau massQ Backgrounds from tau decays to fi- 



^ Corrections for r event feed-across into the background 
estimator are negligible. 
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Table 1 



Summary of 5 ft. and 4/i analyses 



Mode 


57r 


37r27r" 


4h 


Total events 


266 


207 


29 X 10^ 


Events in Fit Region 


36 


19 


17 X 10^ 


Fit Region Purity (%) 


99 


93 


90 


Selection Efficiency (%) 


3.1 


0.4 


2.6 


Typical Mass Resolution (MeV/c^ ) 


15 


25 


20 


Typical Energy Resolution (MeV) 


25 


50 


38 


Uncorrected Upper Limit @ 95 % CL (MeV/c^ ) 


31 


33 


26 
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Figure 1. The number of events in the the 
Fit Region relative to the number in the Con- 
trol Region in the scaled hadronic energy ver- 
sus hadronic mass plane is a function of neutrino 
mass. Kinematically allowed neutrino mass con- 
tours are drawn for neutrino masses of 0, 30, 60 
and 100 MeV/c^ . Note the typical error ellipse 
drawn in the fit region. 



nal states other than 57r and 371271*^ are estimated 
from Monte Carlo studies^, and are expected to 
be small. The hadronic energy scaled to beam en- 
ergy versus hadronic mass distribution is shown 
in Fig ^ for events in the neutrino mass sensitive 
Fit Region. 

4. The Ah Sample 

The scaled energy versus mass distribution of 
the full sample of 4/i decays is shown in Fig ^. 
As noted in table there about 29K events in 
the sample, with 17K events in the fit region. In 



addition to the dominant 37r7r'^ final state, the 
data sample includes other 3ft7r*' resulting from 
the KsTttt'^, Klu, KKsTt'^ and Trcjvr" intermediate 
hadronic states. All these are included in the fit. 
The backgrounds are predominantly from misre- 
constructed tau events in the Fit Region (7%) 
and decays of qq pairs (3%). These background 
shapes are estimated from Monte Carlo. The tau 
background sample normalization is estimated 
from Monte Carlo, and the qq background nor- 
malized to the number of events seen in the data 
above the tau mass, corrected for r feedthroughs. 



5. Spectral Functions 

The invariant mass spectra used in the 5h anal- 
ysis are fit in an independent data sample of tt 
tagged 5h decays. The fit functions are derived 
from that expected from e^e^ — > 47r distribu- 
tions, along with soft pion theorems, as well as 
the shape expected for t GTru, as extrapolated 
from e+e~ 6n distributions |S-10|. The for- 
mer shape is expected to dominate, while the lat- 
ter shape is allowed in the fit as a purely phe- 
nomenological tuning parameter. Since the Qtt 
shape expects events at higher masses than the 
Ait shape, its inclusion naturally leads to a more 
conservative limit. The fit also includes a back- 
ground shape, and excludes the mass range within 
100 MeV/c^ of the tau mass to decrease bias on 
the neutrino mass. The final distributions are 
shown in Fig ^. The 4/i spectral function is dom- 
inated by the 37r7r° final state. This state arises 



2 The CLEO Monte Carlo is based on 



via intermediate decays including a, p, p' , p" or 
an UJ. The spectral function is obtained by the 
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Figure 2. The hadronic mass for the Stt (a) and 
37r27r° (b) event candidates. The Monte Carfo 
estimation of signal shape, normahzed to the data 
sample size in the control region, is displayed as 
the dotted histogram. The solid histogram is the 
background estimator displayed at five times its 
true size for illustration purposes. 

fitting all 27r, Btt and 47r combinations in the cur- 
rent to a sum of such intermediate states, as an 
extension of the model used in Tauola[^: 

+ F--(gL)s.=i,2A/rfe) 

where F''^'^, F'^^ , and are sums of Breit- 
Wigner distributions for the p, p' and p", and 
is a Breit-Wigner distribution for the u. The 
47r combinations used are only those with a mass 
below 1.6 GeV/c^ . Other 3/i7r° decays modes 
not considered in the above are small and fixed 
to their Monte Carlo expected shapes. Although 
this is not a unique model for this decay, all mass 
distribution projections fit very well - see Fig ^. 

6. The Likelihood Functions 

The event likelihood includes both a signal 
shape with a dependence on neutrino mass, and 



Figure 3. The hadronic energy (scaled to beam 
energy) vs hadronic mass for (a) the Stt and the 
(b) 37r27r° event candidates in the fit region. El- 
lipses represent 68% confidence levels. 



a background shape: 

Cevent{M^) = aCsig{data, M„) 
+ {I ~ a)Cbgd{data) 

where a sets the relative background and signal 
size f5 The full likelihood is a product over the 
individual event likelihoods, and, for the 5/i anal- 
ysis, includes a Poisson coefficient P{Nobs, M^i^) 
to account for the expected number of events in 
the Fit Region given the number of events in the 
Control Region, and the neutrino mass. 

The signal likelihood used in previous analy- 
ses relied on convolutions of parametrizations 
of Monte Carlo data. The present analyses forego 
the unnecessary convolution and parametrisation 
steps, and fit directly to Monte Carlo data. Re- 

^ Note that an a priori assumption on the neutrino mass 
allowed range must be made for normalization purposes. 
In this case we chose to assume that Mi, < 100 MeV/c^. 
The result is not strongly dependent on this assumption 
for reasonable a priori distributions. 
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Figure 4. The hadronic energy (scaled to beam 
energy) vs hadronic mass (scaled to r mass ) for 
the Ah sample. The solid lines indicate the Fit 
Region, while the dashed line indicates a larger 
fit region used in systematics studies. 



constructed Monte Carlo events, generated with 
a massless neutrino, are accepted or rejected by 
the data selection criteria. Their generated mass 
and energy are then smeared according to an ana- 
lytic detector smearing function, depending only 
on hadronic mass and energy, and their associ- 
ated propagated tracking and calorimetry errors. 
Furthermore, this smearing function is taken as 
universal for all events in the same sample when 
scaled according to propagated mass and energy 
errors. The probability that a Monte Carlo event 
could be reconstructed as the data event in ques- 
tion is proportional to the smearing. Symboli- 
cally, 

oc J2 G( ^'°*°'^'"^ )>V(Xmc, M,) (1) 

MC 

where G represents the detector smearing func- 
tion, and W is an exactly calculable weight factor 
for neutrino mass hypotheses relative to a mass- 
less neutrino. The smearing function is obtained 
from Monte Carlo r decays and is parametrised 
as the sum of three 2 dimensional Gaussians. 
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Figure 5. The fitted spectral functions compared 
to the data for the 5h samples. The top (bottom) 
plots show the 37r27r° (57r) sample mass fit and re- 
sulting energy projections. The crosses represent 
the 5h data, while the dots represent the fit, while 
the estimated errors in the fit are shown by the 
histograms. Note that the fit was performed on a 
statistically independent sample. 



7. Results 

The final likelihoods are shown in Fig |^ and 
Fig ||. An integral of the likelihood from 
MeV/c^ to 95% of the full area of the likelihood 
gives raw upper limits of 33, 31 and 26 MeV/c^ 
respectively for the 37t27t'^, Stt and Ah samples. 
Combining the 3tt2tt'~' and Stt likelihoods give a 
raw upper limit of 27 MeV/c^ . Including the 
Ah likelihood lowers the uncorrected upper limit 
to 25 MeV/c^ . Systematic errors were care- 
fully examined for the combined 5h studies and 
the Ah study. Charm meson decays to Ktt, K2w 
and K3tt final states in data and Monte Carlo 
show good agreement for the smearing functions, 
and are used to estimate systematic errors for the 
smearing function along the mass axis. Since B 
mesons at CLEO are produced at rest, the mass 
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Figure 6. The spectral function used in the 
4ft. sample fit, for events containing no lu (top), 
those containing an uj (middle), and all events 
(bottom). 



Figure 7. The likelihood function versus neu- 
trino mass from the 5h sample fit for the 871271" 
sample (top), Stt sample (middle) and the com- 
bined 5/1 sample (bottom). 



spectra of hadronic final states of B meson de- 
cays give a direct measure of the smearing along 
the energy axis. Smearing modeling systematics 
contribute 1.5 and 0.4 MeV/c^ to the 5h and 4/i 
studies respectively. The spectral functions are 
all allowed to vary according to the uncertainties 
in their associated fits, contributing 1.9 and 1.2 
MeV/c^ respectively. The CLEO mass and mo- 
mentum reconstruction scales are also considered 
as sources of error (1.5 and 2.3 MeV/c^ respec- 
tively) as was the CLEO energy scale (0.2 and 
3.7 MeV/c^ ) ^ The sum in quadrature of all con- 
sidered systematic errors is 3.1 MeV/c^ for the 
5h study and 5.1 MeV/c^ for the Ah study. Fol- 
lowing the conservative prescription used by the 
LEP experiments |@], these errors are added lin- 



* Differences in nomenclature for systematic errors be- 
tween the studies explain the apparent size inconsisten- 
cies of systematic errors. For instance, modeling of the 
contribution in the smearing is accounted for as a smear- 
ing systematic in the 5h study, and as a momentum scale 
systematic error in the 4/i study. 



M^" MeV/c^ 


25 events 


450 events 





3% 


67% 


50 


« 1% 


< 1% 



Table 2 

Probability of a 95th percentile below 27 MeV/c^ 
for small and large data samples given a true neu- 
trino mass of and 50 MeV/c^ . 



early to the upper limits, resulting in an upper 
limit of 30 MeV/c^ for the 5h sample and a pre- 
liminary limit of 31 MeV/c^ for the 4/i sample, all 
at the 95% confidence level. The combination of 
the 5/1 and Ah studies is still under consideration 
at this writing. The interpretation of the 95th 
percentile of a likelihood as a statement of proba- 
bility about the neutrino mass is dependent upon 
the experimentally reconstructed likelihood being 
representative of an ensemble of likelihoods for 
samples drawn from a distribution with similar 
properties. Thus conclusions drawn from exper- 
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Figure 8. The likelihood function versus neu- 
trino mass from the Ah sample fit. The small bin 
to bin variations result from limited Monte Carlo 
statistics, and are not significant. 



iments with small statistics must be interpreted 
with care especially when very stringent limits are 
quoted. These limits tend to have much lower dis- 
criminatory power than one might naively think. 
Table ^ illustrates this for studies done with the 
CLEO Monte Carlo. For a sample of 25 events, a 
95th percentile limit of 27 MeV/c^ is obtainable 
3% of the time with a massless neutrino but it 
is also obtainable about 1% of the time with a 
50 MeV/c^ neutrino. The discriminatory power 
is of course much better for larger samples. Ex- 
tensive Monte Carlo studies have shown that the 
95th percentile upper limits obtained for the 5h 
samples in this study are not unlikely - in some 
35 % (74 %) of Monte Carlo experiments with 
similar statistics, the 95th percentile was found 
to be larger than that found in the data for the 
5tt ( 3tt2tt'~' ) sample. Thus the interpretation of 
the combined 95th percentile as an upper limit is 
credible. 

In all the fits presented here, the likelihood 
peaks away from MeV/c^ . In no case how- 
ever is the peaking deemed to be significant. In 
the 4h study, the ratio of likelihoods at the peak 
relative to that at MeV/c^ corresponds to a 1.6 
a effect if one assumes that the underlying prob- 
ability densities are Gaussian. Multiple Monte 
Carlo experiments, as well as basic statistical the- 



ory, also show that the likelihood function is not 
always expected to peak at MeV/c^ , even for a 
truely massless neutrino. In addition the conser- 
vative approximations made in the analysis can 
not only result in less stringent limits, but can 
also have the effect of pushing any peak further 
from MeV/c^ . These likelihood curves should 
not therefore be used to infer the existence of a 
massive neutrino. 

8. Conclusions 

Using the world's largest dataset, CLEO sets 
a combined upper limit of 30 MeV/c^ for a com- 
bined study of T ^ 57ri/, r — > 37r27r°i/ decays and 
a preliminary upper limit of 31 MeV/c^ from the 
decay r ^ Shv. These limits use a novel Monte 
Carlo integration technique, and are the first to 
explicitly account for possible backgrounds in the 
likelihood function. 
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